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Abstract 

We propose a two-step procedure to study the order of phase transi- 
tions at finite temperature in electroweak theory and in simplified models 
thereof. 

In a first step a coarse grained free energy is computed by perturbative 
methods. It is obtained in the form of a 3-dimensional perfect lattice action 
by a block spin transformation. It has finite temperature dependent coeffi- 
cients. In this way the UV-problem and the infrared problem is separated 
in a clean way. In the second step the effective 3-dimensional lattice the- 
ory is treated in a nonperturbative way, either by the Feynman-Bogoliubov 
method (solution of a gap equation) , by real space renormalization group 
methods, or by computer simulations. 

In this paper we outline the principles for (/j 4 -theory and scalar electro- 
dynamics. The Balaban-Jaffe block spin transformation for the gauge field 
is used. It is known how to extend this transformation to the nonabelian 
case, but this will not be discussed here. 
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1 Introduction 



The possibility of restoration of spontaneously broken symmetry in the elec- 
troweak theory at high temperature has recently led to a renewed interest in the 
electroweak phase transition JT|. 

Since the seminal and pioneering paper by Kirzhnits and Linde 0, consider- 
able effort has been devoted to describe the rather involved processes occurring 
very close to the critical temperature. The effective potential has shown great 
usefulness. It gives the free energy as a function of the magnetization. There 
are perturbative computations combined with 1/N expansions for the effective 
potential ||, and also numerical results based on Monte Carlo simulations and 
on three dimensional reduced actions 

Nevertheless the perturbative calculations have shown to be plagued with 
problems which are to a large extent due to infrared divergences and are man- 
ifested in the appearance of spurious complex terms in the expansion for the 
effective potential. A great amount of work has been done to obtain meaningful 
expressions very close to the critical surface. Most of these attempts are related to 
the resummation of an infinite number of daisy and superdaisy diagrams which 
have to be taken into account in a consistent perturbative expansion, because 
they are of the same order of magnitude in the coupling constants [||]. 

Another very successful approach was proposed by Buchmiiller et al. || , who 
use an improved perturbation expansion for the effective potential, where the 
dynamically generated plasma masses are already included in the corresponding 
expressions from the beginning. These plasma masses are computed up to one 
loop order by self-consistent gap-equations, and have the effect of damping the 
infrared divergences. In this context, in the absence of linear terms in the 
scalar field of the effective potential was explicitly verified for the abelian Higgs 
model, up to the order e 4 and A 2 , and proved to remain valid at higher orders. 
Such a linear term contains spurious infrared divergent contributions. 

Despite the important and promising contributions to the understanding of 
the nature of the phase transition that have already been made, there are still am- 
biguities to be explained and higher order corrections to be included consistently 
to ensure the survival of the predictions to all orders. 

We propose an alternative variational two-step method to study phase tran- 
sitions at finite temperature in the electroweak theory, and in related models. It 
consists first in the calculation of a coarse grained free energy by perturbative 
methods, obtaining a 3-dimensional perfect lattice action as a function of a block 
spin field. Their coefficients depend explicitly on the temperature . In a sec- 
ond step, the Feynman-Bogoliubov method can be used with this 3-dimensional 
lattice action to obtain the best quadratic approximation to the perfect action, 
and therefore the best split into a free part and an interaction. Using this, the 
constraint effective potential can be computed. To lowest order, the masses are 
obtained in this way as solutions of a gap equation whose formal solution is the 
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sum of superdaisy diagrams for the 3-dimensional lattice theory. They depend 
on temperature. Higher order corrections can be computed in principle. These 
masses are to be inserted into the expression for the constraint effective potential. 

This procedure separates the UV- and IR-problems in a clean way. It allows to 
study the perfect 3-dimensional lattice action directly by numerical simulations. 
One could for instance compute the constraint effective potential numerically. 

The paper is organized as follows. In section 3 the block spin transformation 
for the scalar fields at zero temperature is discussed. The evaluation of the 
perfect scalar action at T=0 is explained and the extension to finite temperature 
is outlined. An explicit computation of the leading terms of the perfect action 
for the A<3> 4 — theory is shown here as well. In section 4 the Feynman-Bogoliubov 
method is applied to the lattice action of the previous section and a comparison 
with other approaches is performed here and in appendix [B]. In section 5 we 
extend the procedure to the Maxwell theory at zero temperature. In particular, 
the averaging operator for abelian gauge fields, the perfect lattice action, the 
gauge fixing and the interpolation kernel are discussed. The block spin of Balaban 
and Jaffe for the abelian gauge field is used ||. In section 6 the procedure is 
extended to scalar electrodynamics. Here the main features of the extension of 
the perfect lattice action to finite temperature are outlined and some quantities 
which are needed for perturbative calculations are exhibited. A short version of 
the results of this paper was presented at LATTICE 94 in Bielefeld 0. 



2 Scholium : Free propagators at finite temper- 
ature 

To prepare for the block spin transformation at finite temperature, we recall first 
the well known relation between free propagators at zero temperature and at 
finite temperature T = 

Because of translation invariance, the standard free propagators in Euclid- 
ian space time depend on a single variable x = (r, t). The zero temperature 
propagator has the form 

v (r,t) = / — — - — exp(-zp ■ r - iut) uj + p + m (l) 
J {Zir) ,) Itt l j 

A finite temperature field theory lives on a Euclidean space time which is 
periodic in time direction with period (3, i.e. on a tube. Propagators admit 
a well known random walk representation. The difference between zero and 
finite temperature stems from the possibility that these random walks may wind 
around the tube several times. Accordingly, the finite temperature propagator 
vt is obtained from vq by periodizing in time. 
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v T (r,t) = 5> (r,t + n/3) (2) 

neZ 



One may use Poisson's resummation formula. 



exp(— iuj{t + n(3)) = exp(-iujt) V] 27T/3 5(u; — cj n ), with cj„ = 2nnP 1 (3) 

As a result , the well known finite temperature propagator is obtained as a sum 
over Matsubara frequencies u n 

— exp(-zp ■ r - iu n t) [uj 2 n + p 2 + m 2 J (4) 



Later on, we will concentrate on the m = case. 

3 Scalar fields only 

To introduce our method we start with a theory with scalar fields only. We 
consider first the T = case and define the blockspin, the effective action, the 
interpolation operator and the fluctuation field propagator. Then we deal with 
the T > case i.e. with the generalizations of the above mentioned entities. As 
an illustration we write down the effective action for T > to leading order in 
the fluctuation field propagator. Finally we add a note on the cancellation of the 
T-dependent UV divergencies. 

3.1 Block spins for scalar fields at zero temperature 

We start from a scalar field theory on the continuum i? 4 of points z. The con- 
tinuum is divided into blocks x. Normally one chooses hypercubes, but here we 
wish to admit different extensions L s of the blocks in space direction and L t in 
time direction. We identify the blocks x with the sites at their centers. In this 
way a lattice A of lattice spacing (L s , L t ) in space and time direction is obtained. 
With a scalar field f{z) one associates a block spin . Following Gawedzki 
and Kupiainien JT6J, we choose them as block averages. 



= CV(x) = av ip(z) (5) 

C is called the averaging operator. It has a kernel C(x,z) which equals the 
properly normalized characteristic function Xx( z ) °f the block x. Because of in- 
variance under simultaneous translations of z and x by lattice vectors in A, its 
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Fourier representation has the form 
C(x,z) = j^—Xxiz) 

-j—^e^{-ip{z-x)-ilz)C{l,p) (6) 

Summation over I = (1, I4) runs over Z4 G 2nL~ l ~ 1 Z and 1 G (2nLj 1 Z) 3 . Integration 
over p is over the first Brillouin zone. 

IIpBvtL; 1 , |p 4 |<7rL t - 1 (7) 

Explicitly 

cu^ = a(z 4 ,P4)a(i,p) 

C t {h,PA) = ~n — 7T~\ sm ' 



L t (p 4 + l 4 ) 2 
3.2 Definition of the perfect action 

Given the continuum action L((p) , the perfect 4-dimensional lattice action La($) 
at zero temperature is defined in Wilson's way 

exp(-L A ($)) =jD<p 5{Cip - $) exp(-L(^)) (9) 
In place of the ^-function, one may use a Gaussian 

6 K (C<p - *) = II (— )" S exp(-^(C^(o;) - $(a;)) 2 ) (10) 



Hasenfratz and Niedermayer |T7| pointed out that optimal locality properties 
of the perfect action L\ are obtained for a preferred finite value of k. 

3.3 Evaluation of the perfect action, interpolation opera- 
tor and the fluctuation field propagator 



To evaluate L\(<&) by perturbation expansion one proceeds as follows One 
splits 

L(<p) = L ((p) + V(<p) (11) 
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If 1 

LoW = 2 J % ( P d iM> = -^(V.^) ( 12 ) 

It is convenient to include mass terms in the interaction V. 

Given the block spin $, one determines a background field if) which minimizes 
Lq subject to the constraint Cif) = $. Because L is quadratic, if> is a linear 
function of $. 

if;(z) = A<S>{z)= [ A(z,x)$(x), with C^l = 1 (13) 



and 

/ =L t L 3 s J2' x = (nL s ,n 4 L t ). (14) 

In coordinate space 

d-U-CV 1 , u = CvC\ v = (-A)- 1 (15) 

Following Gawedzki and Kupiainen [[l(J one splits the field into a low fre- 
quency part, which is determined by the block spin $ , and a high frequency or 
fluctuation field C which has vanishing block average. 

if = A® + C, C( = (16) 

Because of the extremality property of if), the kinetic energy splits 

Lo(<p) = \{A$, -AA*) + \((, -AC) = \{®, u- 1 ®) + i(C, -AC) (17) 

From this one sees that the free massless propagator for the block spin $ equals 
u. Since CA = 1, we have S K (dp — $) = <5 K (CC)- 
The measure splits 

Z^ K (CV-$)exp(-LoM) = ^^ 

The C-dependent factor is proportional to a Gaussian measure with covariance 
(free propagator) T 

D(8 K (C() exp(-i(C, -AC)) = Z-^MC) , r = (-A + k&C)- 1 (19) 



The limit k — > oo can be taken if desired, and results in [If 

r = u - .a^a* = u - ACucW , cr = rc f = o (20) 
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The interpolation kernel A(z, x) of A has a Fourier expansion like C in eq. (|6|) 
with kernel 

A(l,p)=d(k)C*(l,p)u- 1 (p) (21) 

u (jp) = y^ y v(k) I C(l,p) | 2 = / u(x) exp(— ipx) (22) 

We use the abbreviation k = p + I here and throughout. We note that 

u(p) ~ \ for p 2 -> (23) 

because only the / = term in ( |2"2"D is singular at p — 0. 

The fluctuation field propagator r(z, 2; ) is invariant under translations by 
lattice vectors in A. It admits therefore a Fourier expansion of the form 

Summations and integrations are as explained after (Q). All the summations over 
I converge well. For k = 00 , eq. (p0|) yields the explicit expression 

T(l,p,l')=v(k)Si,i'-M,PHp)A*(l',p) (25) 

with k = p + I, and ^4, ii from eqs. fl2"T|) and (fZ2j). 
Since Z is constant, it follows that 

Leff($) = 7i(®> U ^®) + ^//( $ ) + COI1St - ( 26 ) 

exp(-K//($)) = / ^r(0exp(-V*(0) with V*{C) = V (A® + () . (27) 

Thus, V e ff is the f ree energy of a field theory with free propagator T and 
^-dependent coupling constants. It admits a standard perturbation expansion 
(cumulant expansion). We will write down the leading terms for the T > case 
below. All this is well known. 

3.4 Scalar fields at finite temperature 

At finite temperature, we have periodicity in time with period (3. The extension 
L t of blocks in time direction must be chosen commensurate with /3. 

A great simplification results if we chose L t = /3, so that only one block fits 
in time direction. The lattice A becomes a 3-dimensional lattice. 
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Our 3-dimensional fields still have the same dimension as the original 4- 
dimensional ones. This could be remedied by a rescaling by (3^. We prefer not 
to do so because eq.(|51D below looks very natural. Since we do not rescale the 
block spin $ in position space the appropriate 3-dimensional integration measure 
contains a factor (3 and the 3-dimensional 5-functions a corresponding factor 
as follows. 



x=(nL SJ 0). (28) 



In momentum space it is correspondingly. 

To adjust to the finite temperature situation, we need to periodize the T = 
quantities in time in the manner described in section 2. Let us begin with the 
fluctuation field propagator of eq. fl2"DD . We write z — (z, t). 



r T (z,t;z,0 = £r(M + n/3;z,t) (29) 

We insert the Fourier expansion for V . Since I4, l 4 G 2nL^ 1 Z = 2ir(3~ 1 Z, 
we have / 4 t = l±(t + n/3) mod 2n. Only a sum I]n ex P(^P4(^ + n P)) needs to be 
done. This sum was evaluated in section 2. As a result 



r T ( z, t; z', = S 4 / —^3 exp(zk ■ z - ik! ■ z' + il A t - il' 4 t')T(l, p, 0, V) (30) 
l,V P J (27r) 

where k = p + l and k' = p + l'. 

In other words, periodization has the effect of setting p 4 = 0. This reflects the 
fact that the lattice A is only 3-dimensional. 

Similarly one finds the interpolation and averaging kernels and the block prop- 
agator. 



A T (l, P ) = A(l, P ,0), C T (l, P ) = C s (l, p)6 hfi (31) 

The block propagator in momentum space is w(p, 0). It is temperature dependent, 
because the range I4 G 2nL~^ 1 Z of the Z 4 -summation depends on (3 — L t . For the 
same reason, the fluctuation field propagator and the interpolation kernel are 
T-dependent. 

It turns out that the temperature dependence of this block spin propagator 
is a simple factor f3~ l . Therefore we agree to extract this factor, writing 

u T {p) = u(p, 0) = pr Vt (p) (32) 
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with T-independent 'finite temperature' propagator uft (p)- In coordinate space 
it is expressed in terms of the original bare zero temperature propagator 

«ft(x — y) = av av / dt vj>(z — z' , t) = av av / dt v (z — z' , t) (33) 

zG[x] z'G[y] JO zG[x] z'G[y] J-oo 

[x] is the 3-dimensional cube of sidelength L s with central point x. 

3.5 Leading terms in the perturbation expansion of the 
perfect action for (^ 4 -theory at finite temperature 

The perfect lattice action at finite temperature equals 

LeffA®) = f 0, «ft$) + K//,t($) (34) 
VeffA®) = - In (/ d^r T (C) exp {-V %T (C))) (35) 



with 



14,t(C) = V (At® + C) (36) 



The effective interaction (including mass terms) is temperature dependent be- 
cause the fluctuation field propagator IV and the interpolation kernel At are 
both temperature dependent. This temperature dependence is weak and disap- 
pears in a zeroth order local approximation. 
Let us consider y9 4 -theory with bare mass m . 

L ((p) = g / {dtfdpip) dz (37) 



V(ip) = J ^-m^ 2 + ^jV^ 4 ^ dz + wave function renormalization term (38) 
Inserting the split fll6|) of the field (p , we obtain 



7.(0 = U d ($) + E [dz -g n (<S>, z) C(z) n (39) 

n=l J U - 



with 
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9i — 


^{zf + m 2 r^ (z) 


92 = 


^(zf + ml 


93 = 


A *{z) 


9a = 


A 


u d ($) = 


/ O + S 


9(z) = 


/ .A T (z,a;) $(z) 



(40) 



V e ff t T can be calculated in a loop expansion. When the starting point is a lattice 
theory, one can use Mayer expansions instead. They are convergent for weak 
coupling, and the asymptotic expansion of individual terms in powers of the bare 



coupling constant contains infinite sets of diagrams |Le 
We write the perturbative expansion of V e ff t T as 



K//,t($) = M$) 



(41) 



N>1 
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where V^f) T scales as 7^ when the fluctuation field propagator IV 

Let us compute v}j\ T . Indicating factors \1> by dotted external lines and a 
fluctuation field propagator by a solid line we have 



Figure [3] effective potential to first order in the fluctuation propagator 



V. 



(i) 



eff,T 



1 

2! 



A 



dz(-^(z) 2 + m 2 )T T (z,z) 



1 
2! 



dz 1 dz 2 (^(z 1 ) 3 + m 2 ^(z 1 ))T T (z 1 ,z 2 )(^(z 2 ) 3 



3! 



3! 



m 2 *(z 2 )X42) 



The fluctuation field propagator T T (z,z), decays exponentially with distance 



z — z 



z — z 



with decay length one block lattice spacing. That is, it decays with 

with decay length L s . Similarly, At{z, x) decays exponentially in \z — x\ 
with decay length L s . As a result, each term in V e ff,T is local in $ modulo 
exponential tails with decay length of one lattice spacing. 
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By a process of partial integration (or rather summation) one can exhibit each 
term in V e jf,T as a sum of local terms and small remainders which represent the 
exponential tails (see appendix |A[). The terms with coefficients of dimensions up 
to (mass)" 2 are 



V, 



eff,T 



($) 



l |im 2 $ 2 + 5 Z (V,$) 2 + + + 7 (V,$) 2 $ 2 + • • •} 



(43) 



Since u F \~{p) ~ p 2 as p — > 0, we may write the expansion equally well in the form 
K//,t($) = 

R r r 

$(x)w^(x, y)$(y) 25 z + 7$(x) 2 



2 4! 6! 



+ 



(44) 



All the coefficients are T-dependent and finite. 

For computer simulations it is appropriate to consider V e jf,T as a function of 
$(x). For analytical computations it is most convenient to regard it as a function 
of*(Y). 



3.6 Note on the cancellation of temperature dependent 
UV-divergent diagrams 

The UV-convergence of a quantum field theory concerns its local behaviour and 
is therefore not temperature dependent. If the proper choice of counter terms 
makes the theory finite at zero temperature, then also at finite temperature. The 
cancellation occurs order by order in perturbation theory. This is well known. 

Problems can occur when one sums selected classes of diagrams. We do not 
propose to do so when deriving the perfect action, but stay strictly within the 
realm of standard perturbation theory. So there can be no problem. 

It is nevertheless appropriate to point out that there do exist individual dia- 
grams with temperature dependent UV-divergent pieces. They cancel. We give 
an example. 

Let us split the (fluctuation field) propagator V into a static part (heavy line) 
r siai and a non static part T ns (wavy line). The non static part represents random 
walks which wind several times around the tube; it is not singular at coinciding 
arguments : T ns (z, z ) < oo . But it is not zero and is T-dependent. As a result, 
the diagrams 





Figure El temperature dependent divergent diagrams 
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are both logarithmically UV-divergent, with a temperature dependent coefficient 
T ns (z,z f ) . They must cancel. SX is the logarithmically divergent coupling con- 
stant counter term to second order in A . 



4 The Feynman-Bogoliubov Method 

To obtain information on the nature of a phase transition at finite temperature, 
one wants to compute the effective potential, i.e. the free energy as a function of 
the magnetization. 



M = J dzLp{z) = J $(ac) (45) 

Alternatively, one may compute the constraint effective potential which gives 
the probability distribution of M. 

To do so we should apply nonperturbative methods to solve the lattice theory. 
One of these methods consists in the solution of gap equations. It is very old and 
known as the Feynman-Bogoliubov method [JT2]] . 

The gap equations have "perturbative" solutions which come from summa- 
tions of superdaisy diagrams. In principle they might also have other solutions. 
We wish to examine this method in order to see what will be the effect of terms 
like (V M $) 2 $ 2 etc. in the perfect lattice action. 

Let us make it clear that it is not sufficient to find the solutions of these gap 
equations. In order to justify the perturbative calculation of the perfect action 
it will be necessary to investigate also the stability properties of the solutions 
of the gap equations against small perturbations of the lattice action. We will 
comment on this, but a thorough treatment of this question is beyond the scope 
of this paper. 

Given an action S(<&) of some Euclidean field theory, one seeks an optimal 
quadratic approximation Sf ree ($>) around which to expand 

S f reem = \l [ ${x)J{x,y)<S>{y) (46) 

A J x J y 



By Peierls inequality |I2|, the partition functions obey the inequality 



\n{Z) > ln{Z free ) - (S-S free ) free (47) 

for any choice of the Sf ree . Herein, <>/ ree is the expectation value in the theory 
with action Sf ree . The right hand side is the first order approximation to the 
left hand side in the perturbative expansion around Sf ree . The optimal choice of 
Sf ree is that which makes the right hand side of ( (4"T| ) maximal. In other words, it 
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makes the first order approximation as good as possible. It is asserted that there 
exists always a unique optimal choice of Sf ree . It is not asserted that the optimal 
choice is necessarily a good one. The optimal J is determined by the extreme 
value condition 



5 2 S \ 
,<5$ (x) 5<5>(y)/ 



= J(x,y) (48) 

free 



This is equivalent to the condition that the right hand side of ( [47|) is maximal, 
viz 



jj [HZfree) ~ (S - % ee ) /ree ] = (49) 

When applied to the standard $ 4 action, this produces the gap equation whose 
perturbative solution is the sum of superdaisy diagrams (see below). 

Let us consider the gap equation which results from the more complicated 
lattice action. 



5($) = |($,^$) + K//,t($) 



(50) 



with V e ff,T from equation (§4] 
We obtain 



5 2 S($) 



5$ (x) 5$(y) 



pup^(x, y) \l + 25 z + 37$(x) 



+ P5(x-y)3'y$(x)(u]?j r $)(x) 

+ 5(x-y)L 2 + ^(xr + ^(xr 



The expectation value in the theory with action Sf ree = J |$J$ is 

= (3u- F l T {x,y) [l + 25, + 3 7 J- 1 (0)" 



' 6 2 S($) \ 
^(x)5$(y)/ free 



+ p5(x-y)3 1 (up 1 T J- 1 )(0) 
+ 5(x-y)(m 2 + ^-J-\0) + 



2 2 2! 



(51) 



J-\0) 2 ) (52) 



where J l (0) = J 1 ( independent of x by translation invariance, assuming 

that we seek a translation invariant solution. Other solutions could be of interest. 
The gap equation fl4"8]) can be solved by the Ansatz 
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J(x,y) = Au F r(x,y) + BS(x - y) 



(53) 



Inserting the Ansatz results in two transcendental equations for A, B whose 
solutions depend on the coefficients (3, 5 Z , 7, Aj. 

Basically, the inclusion of the (V M $) 2 $ 2 -term results in a system of two 
equations for mass and wave functions renormalization. In standard $ 4 -theory 
there is only one equation for the mass. 



4.1 Comparison with other approaches 

The standard method of dealing with the finite temperature phase transition of 
$ 4 -theory is described by Kapusta [|TJ. It consists in the summation of ring or 
daisy diagrams which contribute to the self energy part II of the physical finite 
temperature propagator. 

It would seem natural to improve this by solving the gap equations instead. 
This would replace the sum over daisy diagrams by a sum over superdaisy dia- 
grams. The temperature dependent mass m 2 (T) would be determined in a self 
consistent way. But there is a problem with this: There arises a temperature 
dependent logarithmically divergent contribution to the self energy which can- 
not be cancelled by a temperature-independent bare mass term. The approach 
outlined in this paper avoids this problem, because UV-divergences are cancelled 
before one derives the gap equation. We see here the advantage of making a clean 
separation between UV- and IR-problems. 

Let us explicitely show how the problem arises in a conventional superdaisy 
approach without a lattice as an UV-cutoff. 

We use the standard <3> 4 -action in place of (|50D and make the usual Ansatz 
for the solution of the gap equation. 



J(x,y) = -A5(x-y) + m 2 (T)5(x-y) (54) 

We obtain the gap equation (|4"8| ) as before with 5 Z = 7 = \q = 0, u~ x = — A. It 
is solved if 



m 2 (T)=ml + ^J-\0) (55) 



JT\0) =TE/^R + P 2 + m\T)[ l (56) 
The integral defining J _1 (0) is evaluated in the standard fashion, resulting in 
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j-\o) = n mc + n 



mat 



(57) 



n, 



p 2 + m 2 (T) 



16tt 2 



A 2 — m 2 (T) In 



A 2 



m 2 (T) 



(5? 



n 



mat 



d 3 p 1 

(2tt) 3 w 



[exp(/3u;) - 1] 



-i 



(p 2 + m 2 (T))2 



(59) 



The integrals are UV-divergent, therefore a regulator A is introduced. It is sup- 
posed to be taken to infinity after cancellation against counterterms. We see that 
without a lattice a temperature dependent counterterm is needed to cancel the 
logarithmic UV- divergence in fl58|). This is not acceptable. 

Conclusion: The Feynman-Bogoliubov method is designed to deal with in- 
frared problems. It should only be applied after a correct cancellation of UV- 
divergences has been achieved through the computation of an effective field the- 
ory, e.g. on a lattice, or with another low UV-cutoff. 

We wish to compare our proposal with another approach. It was proposed 
to derive first a 3-dimensional continuum theory. When no UV-regularization 
is introduced, then this 3-dimensional theory is superrenormalizable but not fi- 
nite. It has linear self energy divergences. These divergences are cancelled by 
divergent terms in the 3-dimensional action. Again the problem remains that 
the UV-divergences have not been completely cancelled before nonperturbative 
methods are applied to deal with infrared aspects. In addition, the 3-dimensional 
continuum theory is complicated because of nonlocalities with decay length ^- 
which appear as soon as one goes beyond the leading order in the perturbative 
expansion of the 3-dimensional action (see appendix H). The use of a momentum 
expansion to deal with these nonlocalities would aggravate the UV-problems. 



The 'average action' approach of Wetterich [|13|] is in the same spirit as ours. 
The difference is that we use a lattice. The lattice has the advantage that one can 
do computer simulations. Wetterichs work suggests moreover, that one should 
perform simplifications only on the 1-particle irreducible parts of the perfect 
action (i.e. those not held together by a single T-propagator). It has been known 
for a long time that the 1-particle reducible second term ^IY 1 !/ 3 in eq. (|42 ) is 
important inspite of looking irrelevant fll6|. 



5 Perfect lattice action for Maxwell theory 

We wish to extend the consideration of the previous sections to scalar electrody- 
namics. As a preparation for this we consider the free abelian gauge field. 

In this section we recall the Balaban-Jaffe block spin transformation for the 
free abelian gauge field at zero temperature. 
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5.1 Averaging operators for abelian gauge fields 



Given a vector potential a(z) = a^(z)dz^ on the continuum, we define a block 
spin A on the block lattice. A lives on links b of the block lattice. We use the 
alternative notation 



A[6] = A,(x) (60) 

when b is the link emanating from x in /i-direction. 

We will not distinguish in notation between the averaging operator C for 
gauge fields, and for scalars (gauge transformations). Which one is meant will be 
clear from the context. Thus A = Ca. 

Given a link b from x in /i-direction, and a point z G x, let C z ^ be the straight 
path of length one block lattice spacing in /i-direction starting from z. If fi = 4, 
C z ^ connects z with z + L t e^ and if // ^ 4 it connects z with z + L s e^ (e M = unit 
vector in /i-direction). The blocking procedure is defined by 



A [b] = av a [C Z J (61) 



a[C 2 J=/ dw u a u (w) (62) 

This blocking procedure is covariant under gauge transformations in the following 
sense. If 



a^z) -> a^z) - d^\(z) = a*(z) (63) 

then 

A^x) -> A^x) - V^Aix) (64) 



A = C A (65) 

i.e. 

A(x) = av X(z) (66) 

We use the notation a x for the gauge transform of a , etc. In this notation, the 
covariance property reads 
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Ca x = (Caf x (67) 
We will need the Fourier transform of the C-kernel. Let us write 

A v( x ) = J C fl „(x,z)a„(z) (68) 

The kernel will have a Fourier expansion just like the scalar kernel. Explicitly 

k L 2 k L 

C l u,(l,p) = S lu ,exp{-i- ! Y^) sin^ 1 C(l,p) (69) 

where C(l,p) is the scalar kernel (|J), and no sum over v is implied. L v = L t if 
v = 4 and L v = L s otherwise; k u = p u + l u as usual. 

5.2 Perfect lattice action for the free electromagnetic field 

Using the symbol d for the exterior derivative, da — | (d^a u — d u a^) dz^ A dz u , 
the Maxwell action of the electromagnetic field can be written as 

S M (a) = - < da, da > (70) 

The perfect lattice action associated with this is defined by a formula analo- 
gous to eq. Q27D for the scalar field. There is one difference, however. In order to 
give meaning to the functional integral, some amount of gauge fixing is necessary. 

We wish to obtain a gauge covariant perfect lattice action. Therefore we wish 
to retain the freedom of gauge transformations on the lattice, i.e. one gauge 
degree of freedom per block. Therefore global gauge fixing is not appropriate. 
Instead, gauge fixing is only used locally within each block. 

Consider gauge transformations 



<V — > — A with C X = (71) 

By eq. QS7D , they leave the block gauge field invariant. These are the gauge 
degrees of freedom which will be eliminated by fixing the gauge within a block. 
There remains one gauge degree of freedom per block A(x) = C X(x) which is 
not affected by the fixing. It extends to a global gauge transformation per block, 
X(z) = A(x) for z G x. 

To begin with, a block axial or radial gauge is used. For the purpose of 
perturbative calculations, a transformation to a block Landau gauge is carried 
out later. 
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Suppose for a moment that we start from a theory on a fine lattice A/j ne 
instead of the continuum, a will then live on links / of Af ine . We describe the 
block axial gauge for this situation first. The formal continuum limit will be 
obvious. 

For every block y G A , a maximal tree Ax(y) is selected, for instance the 
comb within each block shown in the figure below. 



Figure |3] maximal trees within 2- dimensional blocks 

Ax(y) consist of links I G Af ine . In the block axial gauge a [I] =0 for all 
I G Ax(y), ye A. 

One defines appropriate ^-functions which are inserted in the functional inte- 
gral 

<W«) = II *(<*[*]) ( 7 2) 

leAx(y) 

5ax{o>) = n $Ax(y){a) (73) 
The perfect lattice action S e ff(A) is defined up to an additive constant by 

exp(-S e //(A)) = / Da5 Ax (a)5(Ca- A) exp(-S M (a)) (74) 

Instead of the block axial gauge, the block radial gauge can be used. In the 
continuum limit it is 

4 

^2(z - x),, a^z) = for z G x . (75) 

5.2.1 Transformation to the block Landau or a-gauge 

For perturbative calculations, the block axial gauge is inappropriate because of 
the bad ultraviolet behaviour of the fluctuation field propagator in this gauge. 
Therefore the Faddeev-Popov method is used to convert to a block Landau or 
a-gauge. 

Sm is invariant under arbitrary gauge transformations while 5 (C a — A) is 
only invariant under those gauge transformations a^a x = a — d\ which obey 
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CA = 



(76) 



C is the averaging of a scalar function over blocks. We wish to eliminate the 
freedom of performing gauge transformations which obey ( [76D by gauge fixing. 

Recalling that d stands for the exterior derivative of a p-form, d* is the 
coderivative. For example d^a = <9 M a M and c^A = if A is a 0-form. In this 
language the Laplacian is given by dd* + d^d. 

To convert to the block a-gauge, one multiplies and divides expression (ff^) 

by 

V(a) = J DX5 (CA) exp {-— < tfa\ >| (77) 

After interchanging the order of integration and performing a gauge transforma- 
tion a —> a~ x in the a-integral 

r 1 1 
exp(— S e ff(A)) = Da——5(Ca-A)exp(-S M (a)-7^-<d^a,d' ( a>) 

D\5(C\)5 Ax (a~ x ) (78) 

The second factor is 1 because A is fixed up to a global transformation by the 
condition 5(a _A [/]) for every link / £ Ax(y) in the maximal tree of block y, and 
the remaining global transformation is absorbed by the 5-function 5(CX). 
It remains to evaluate 



V(a) = J DX5 (CA) exp {-— < <9 f a\ <9V >| (79) 

= J DX5(CX) exp < (A- 1 ^^ - A), A t A(A- 1 9 t a - A) >| (80) 

Thinking of the 5-function as a limit of a Gaussian, one sees that this is a Gaussian 
integral. But since the 5-function restricts the range of the integration the result 
is not just a constant. To show this we introduce the projector R onto those 
functions A which satisfy the constraint CA = 0. R can be chosen orthogonal with 
respect to the scalar product < ■ , A^A ■ >. The integral is split into two parts, 
the integration over (1— R)X is restricted by the <5-function to (1— R)X = whereas 
the integration over RX has no further restriction. Since R is an orthogonal 
projector there exists no mixed term and the integral is a pure Gaussian. 
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V(a) = const exp {- — < d*a, A(l - R)A~ 1 d^a >} (81) 
Inserting this into eq.(|78|) we obtain the final result 



exp(— S e ff(A)) = J Da 5 (C a — A) exp (S M (a) - — < d ] a,ARA- l d^a >)(82) 

The projector R ensures that only the gauge degrees of freedom A with C A = 
- i.e. vanishing block averages of A - are fixed. There remains the freedom of 
performing global gauge transformations on each block. This freedom is reflected 
in the fact that exp(— S e ff(A)) is invariant under gauge transformations A — > 
A A = A — d A, where A is a scalar function on the lattice. 

5.2.2 Explicit formula for R 

We recall from section 3 that the fluctuation field ( associated with a scalar field 
(ft satisfies C ( = and can be obtained by applying a projector 



C=(l-AC)<p (83) 

A C is an orthogonal projector with respect to the scalar product which is fur- 
nished by the kinetic term for (p. 

Here we look for an orthogonal projector with respect to the scalar product 
employing the kinetic term for A which is A 2 in place of —A . Therefore we may 
write 



R=l- A {x) C 



*4) 



where C is the same averaging kernel for scalars as before with its Fourier trans- 
form given by eq.(|6]), while A^ is chosen to satisfy 



v-^AW = C t u -i. 



u 



CvC ] : 



A 



15) 



in complete analogy to eq.(|15D. 
The Fourier transform is 



A w (l, P ) =v(k)C*(l, P ) ^(p)- 1 



(86) 



with 



k — l+ p, 



m = h 



and 



r.A 



U 



(87) 
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5.3 Interpolation kernel A and evaluation of the perfect 
lattice Maxwell action 

The perfect lattice action is determined in essentially the same way as for the 
scalar field. 

Given the block gauge field A, one seeks that gauge field a which minimizes 
Sat (a) +^ < d'a, ARA^ 1 d^a > subject to the constraint C a = A. Because the 
action is quadratic, a is a linear function of A 

a = AA i.e. a^z) = f A b) A [6] (88) 

J links 6eA 

Balaban and Jaffe show that a - and therefore A - are independent of a and 
satisfy the block Landau gauge condition 

Tld j a = i.e. Kd^ A =0 (89) 

where 1Z is a shorthand notation for Ai?A _1 . The proof is based on constructing 
a gauge transformation A with C A = which ensures (p9|). The equation 

A\ = TZd ] a (90) 

has a unique solution in the subspace C A = because 1Z is the projector onto 
those configurations A A which satisfy C A = 0. 

The kernel A^z, b) has exponential decay in the distance of z from the blocks 
(x,y)=b !. 

The effective action is obtained by shifting the field a = A A + a where a is 
called the fluctuation field. The integral over a produces merely a constant and 
one obtains 

exp(-S e// (A)) = exp (~ < &4A, 8AA > -— < d^AA, TZd^AA (91) 

= exp (~2< dAA > dAA >) ( 92 ) 

independent of a because of eq. (|S9D. Clearly S e ff(A) is a bilinear function 
of A. One may therefore define an operator Ai acting on block gauge field 
configurations A, such that 

S eff (A) = ^<A,A 1 A> (93) 

Explicitly 

A l = A ] d ] dA (94) 
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5.4 Fluctuation field propagator 

The probability distribution of the fluctuation field is given up to a normalization 
factor by (« = oo) 

6 K (Ca) exp (— - < da ,da > -— < d ] a,TZd^a > X (95) 

Thinking again of the 5-function as a limit of a Gaussian S K 



8 K {Ca') = (^) 2 exp < Ca', Co' >) (96) 

we see that the probability distribution is given by a Gaussian measure with 
covariance 



T a = lim r Kia (97) 



C is again the Balaban-Jaffe averaging kernel for the gauge field [JTT . 

The block Landau gauge is a = . There exists a formula which expresses 
the fluctuation field propagator T in the block Landau gauge in terms of the 
propagator r KjQ for finite k (e.g. k — 1 ) and arbitrary a ||. 

r = G--Gdftd t G-GC t (CG'C ,t r 1 CG (99) 
where G = r l Q ,. 

The propagators G and r decay exponentially in coordinate space with decay 
length one block lattice spacing. When the original theory lives on a lattice, 
the proofs of reference |J apply. One can also convince oneself of the fact by 
examining the behaviour in momentum space near zero momentum. The explicit 
formulae given below can be used for that. 

The interpolation kernel A (which is independent of a as we know) can be 
expressed in terms of G = Ti a as well 



A =GC ] {CGC^) 1 (100) 

Therefore, the block Landau gauge fluctuation field propagator r can also 
be written as 
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r = G - —GdTZd j G - AC G C ] A ] 

a 



(101) 



This resembles the formula for the scalar case. In place of the propagator for the 
full theory, which does not exist because we have not fixed the gauge completely, 
there appears the auxiliary quantity G = r l a . The expression for T is valid for 
any choice of a in this auxiliary propagator. 

5.5 Perfect Lattice action for the Maxwell field at Tem- 
perature T > 

To go from temperature T = to finite T one must periodize the propagators 
and interpolation kernels in time as in the scalar case. 

D = (-Ax)" 1 (102) 

with Fourier transform D(p). Then the free lattice photon propagator in the 
3-dimensional lattice theory which is obtained from the finite temperature field 
theory is 

D T (p) = D(p,0) (103) 
Similarly, the interpolation kernel 

Ar(Z,p) = .A(Z,p,0) (104) 

and the fluctuation field propagator 

f OT (/,p,0 = f (/,p,0,0. (105) 

As in the scalar case , the discrete variables /, t remain 4-dimensional. Think 
of them as referring to the original 4-dimensional theory. 

6 Scalar electrodynamics 

Now we are prepared to deal with scalar electrodynamics. The definitions for the 
gauge fields can be taken literally from the preceding section. Since the Higgs 
field is not invariant under a gauge transformation we cannot simply use the 
results of section 3. So it is the Higgs field we are concerned with in this section. 
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6.1 Block spin transformation for the Higgs field at zero 
temperature 

Given a Higgs field <ft(z), we wish to define a block Higgs field <&(x) on the lattice 
A in a gauge covariant way. The blocking procedure for the <^ 4 -theory cannot be 
used as it stands because it is not gauge covariant. In order to maintain gauge 
covariance, one must use averaging kernels depending on the gauge field a. 
We denote the averaging operator for the Higgs field by C H (a) 

<S>=C H {a)<j) (106) 



i.e. = f dz C H (a\x,z) <f>(z) (107) 

The adjoint kernel C H ^(a\z,x) = C (a\x,z) 

In the ^-theory we used an averaging kernel which was constant on blocks 
and vanished outside. A constant is the lowest eigenvector of the Laplacian with 
Neumann boundary conditions on block boundaries. 

The natural generalization to the gauge covariant situation is as follows Q. Let 
— A a = (<9 — a)^ (<9 — a) denote the covariant Laplacian and — A^' x the covariant 
Laplacian with Neumann boundary conditions on the block boundary of x. We 
demand that 

- A*' x C H \a\z, x) = e {a\x) C H \a\z, x) (108) 

and 

C H \a\z,x) = for z^x (109) 

where Eo(a\x) is the lowest eigenvalue of —A^' x . In addition we impose the nor- 
malization condition 

C H C H ^ = 1 (110) 
This leaves the freedom of multiplying C H with an a-dependent phase factor 

C H \a\z,x) — >C H \a\z,x) rj(a\x) (111) 
1 Here we deviate from the work of Balaban, Jaffe and Imbrie 
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It follows from the gauge covariance of the eigenvalue problem ( |108| ) that 
under a gauge transformation 



C H (a x )(f) X = (C H (a)(J)) A (112) 

with some gauge transformation A on the lattice which depends on A and on 
the choice of conventions to fix C H uniquely. The freedom in the choice of an 
a-dependent phase factor may be exploited to demand that 



A(x) = C \{x) (113) 



where C is the scalar averaging kernel introduced in subsection |3.1| . Eq. ( |112j ) 



parallels the gauge covariance property of the blocking procedure for the gauge 
field. 

One can compute the averaging kernel C H (a) as a solution of the eigenvalue 
equation ( |108| ) by standard quantum mechanical perturbation theory. A proto- 
type of such a computation is found in reference |14] . 



6.2 The perfect action of scalar electrodynamics: Defini- 
tion 

The perfect action is defined by 

exp (-&//($, A)) = J Da J ' 5 Ax (a) 5 (C a - A) 5 K (C H (a)<P - <$> 



exp (-S M {a) - \ (</>, -A a 0) - V (0)) (114) 



where 



v{<p) 



dz 



1 

-m, 



o 



A 
4! 



(115) 



To evaluate it perturbatively, we proceed similarly as in y) 4 -theory and pure 
abelian gauge theory. 

We prefer to use a Gaussian 8 K {C H {a)(j) — $^ in place of a ^-function. 



6.3 Interpolation kernel for the Higgs field and Higgs fluc- 
tuation propagator 

The interpolation kernel A H for the Higgs field will also depend on the gauge field 
a because the averaging kernel does. In order to split the kinetic term for the 
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Higgs field, we impose the usual demand that ^ = A H $ minimizes — A a \I/) 
subject to the constraint that C H (a)^> = $. This happens if 

A a A H (a\z,x) = ( C H \a\ z,y)A eff (a\ y,x) (116) 

JyeA 

for some choice of A e ff, and if 

C H (a)A H {a) = l (117) 
The last condition implies that 

(A e// (a))" 1 = C H (a) K l C H \a) (118) 

similarly as before (where we used the notation u~ l in place of A e jj). 
If we make a shift of the Higgs field 

( f ) = A H <5> + ( (119) 

The fluctuation field propagator for the Higgs field is 

T H (z, w) = (-A a + k C H ^C H Y l {z, w) (120) 

In the limit k — - > oo it becomes 

T H = v H - A H C H v H C H ^A H ^ (121) 

where v H = — A^ 1 is the full gauge covariant free massless propagator for the 
Higgs field. 

The interpolation operator and the gauge covariant fluctuation field propaga- 
tor can also be computed by standard quantum mechanical perturbation theory. 

6.4 The perfect action of scalar electrodynamics: Rep- 
resentation as a functional integral in block Landau 
gauge 

We use again the Faddeev- Popov trick to convert from block axial gauge to block 
Landau or a-gauge. This is done in exactly the same way as for the pure abelian 
gauge field theory. 

The essential point which makes this possible is that both ^-functions 

S(Ca- A)8 K (C H (a)(f)-<S>) (122) 
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are invariant under gauge transformations a — > a x , — > A which obey the 
constraint C\ = 0. Proceeding as before, we obtain 

exp (S e ff($, A)) = j Da J D(j)5(Ca- A) 5 K (C H (a)<j)-<5>) 

exp (-S M (a) - 5 5/ (a) - ± (0, -A o 0) - V(0))(123) 



where 



S gf (a) = ±-(tfa,Ktf a) (124) 



2a 

with the same 1Z as for the pure abelian gauge field. 
Now we shift the fields as 



= A H $ + C 

a = AA + a' (125) 



The expression for the perfect action becomes 



exp(-S e// (ct>,A)) = J Da 1 j D( 5(Ca>)5 K (C H (a)() 



exp 



- S M ,e//(A) - S af (AA) - S M (a>) - S gf (a') 

1 

2 



~\ -A e //(a) $ ) - i (C, -A a C) - V (A H (a)<S> + () 1(126) 



where a = AA + a' . 

We separate the terms of zeroth order in a and £. 

A e// (a) = A e// (.AA) + E (a, A) = A e// (A) + £(a', A) (127) 
V(A H (a)$ + C) = C/d (A, $) + W(A, a', $, C) (128) 

where 

C/ d (A,$) = \/(^ H (^lA)$) (129) 
The perfect action becomes 
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S e// (<t>, A) = S M , eff (A) + S gf (AA) + U d (A, $) - - ($, -A e// (A)$) 

A) (130) 



To zeroth order in the fluctuation field propagators S e jf(&,A) = after 
subtracting a constant. 

6.5 Scalar electrodynamics at finite temperature 

We wish to adopt the considerations of the preceding section to finite temperature 



We introduce a block lattice A as for the 4 -theory and the Maxwell theory. 
The side length of the blocks in time direction is (3. As a result, the block lattice 
is 3-dimensional. 

The blocking of the gauge field at finite temperature was discussed before. 

For the Higgs field, some new features appear, because the averaging and 
interpolation kernels are gauge field dependent. For this reason, a covariant mo- 
mentum space description does not exist, and we cannot apply exactly the same 
periodization procedure as for a purely scalar theory. Nevertheless, the transition 
to finite temperature by periodization in time is straightforward. The formulae 
of the preceding section remain valid as they stand when they are properly in- 
terpreted. 

First of all, we must regard the blocks as coming equipped with periodic 
boundary conditions in time direction. The Neumann boundary conditions ap- 
ply only to the boundaries of the block which remain, after periodic boundary 
conditions in time direction are imposed. 

Similarly, the defining equation for the interpolation kernel A? remains the 
same as for T = , except that periodic boundary conditions in time direction 
are imposed, and the previously defined averaging kernel is to be used. This 
equation together with the normalization condition C^A^ = 1 also defines the 
propagator 



with 




(131) 



T. 



= -A"/, 



eff,^ 



(132) 
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for the blocked Higgs field. The fluctuation propagator is obtained from these 
quantities as before and will also be temperature dependent. All these quantities 
depend on the gauge field. This implies a hidden T-dependence also for the block 
Higgs field, because the boundary conditions for the gauge field are T-dependent. 

The effective action is defined as before, with substituted for C H , etc.. 
It lives on a 3-dimensional lattice A . Because of the anisotropy, the dependence 
of the resulting action on A and on A,, i = 1,2,3 is different. The A field 
behaves like an extra scalar field in the effective 3-dimensional theory. 
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A Partial integration of nonlocal terms in the 
effective action to obtain a sum of local and 
irrelevant terms 

Consider a term which is quadratic in the field such as 

J 2 ($)=/7 p 2 {x 2 - X!)$(x 2 )$(xi) (133) 

J Jxi,x 2 

We are interested in situations where p 2 (x) falls off exponentially with decay 
length one lattice spacing. In this case the sum can be rewritten in the form 

J 2 ($) = p 2 ®(xi) 2 + <V / V^xi) V^$(xi) + irrelevant term (134) 

J XI J Xl 



where the irrelevant term is of the form 

1p.up, 



■oil ^^uHx l )p' 2 (x 2 -x l )V p V a ^2) (135) 

J JXi,X2 

p' 2 also decays exponentially with distance x 2 — X\. The coefficients are 

/i 2 = / 9%{x) (136) 



V = ~\ I %nXvP2{x) (137) 

Z J X 

Because of the exponential falloff of p 2 , its Fourier transform p 2 {p) is holomorphic 
in a strip. Because of the presence of a lattice, it is a periodic and even function 
of p. Therefore 

h{p) = + z iw sinfV* sin Pf + Tm^po- smpf, sinp u sinp p smp a p' 2 (p) (138) 
where p 2 {p) is also holomorphic, periodic and even. 

p 2 = p 2 (0) (139) 
V = l^^-f^o (140) 

Inserting back one gets eq. ( |131| ). 

B Dimensional reduction to a 3-dimensional the- 
ory in continuous space 

The method outlined in this paper can also be used to compute the dimensionally 
reduced 3-dimensional theory in continuous space by perturbation theory. Its 
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action is local to zeroth order, but will then start to develop nonlocal terms of 
the form 

|p(x 1 ,...,x n )$(x 1 )---$(x n ) (141) 

where p decays exponentially with the length of the shortest tree on vertices 
xi, . . . , x n with decay length 

We outline briefly how this is seen. The method had been used before in the 



work of M. NieBen [19]. He used it to make d-dimensional quantum statistical 
systems more palatable to a computer by discretizing time in the d+l-dimensional 
functional integral representation of the system. 

In our case, the 3-dimensional fields are time averages of the 4-dimensional 
ones. 

$(r) = /T 1 J dtcp(r,t) = J dt'dr'C(r,r',t')(p(r',t') 
C(T,T , ,f)=/3- 1 5(T-i f ) (142) 
The block propagator u(r — r') is therefore given by 

u(y) = (3- l u FT (r) = dtv T {r,t) = /T 1 / dtv (r,t) (143) 

JO J-oo 

it depends on j3 only through the overall factor (note that we did not rescale the 
fields yet). 

The interpolation operator is defined as usual, A = vtC^u^ 1 . Because C is 
time independent, A comes out trivial. 



^(r',t',r) =/3" 1 5(r-r') (144) 

As a result the fluctuation propagator can also be simplified. It is translation 
invariant both in r and t, and 

T T (r, t) = Mr, t) - /T 1 f dt'v T (r, f) (145) 

JO 

The second term can be reexpressed using eq.( |143| ). It is seen that the Fourier- 
transform of IV is the same as of Vt, except for the absence of the zero Matsubara 
frequency mode. Thus 



u n = 2Tmf3~ 1 (146) 
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The decay in r can be read off the singularities of the integrand on the imaginary 
p-axis (p = |p|). The closest singularity is at p = +i2Tr/3~ 1 . Therefore there is 
exponential decay with decay length 

Tt appears as free propagator in the perturbative expansion of the 3-dimensional 
action. Its decay properties will govern the nonlocalities of the resulting action 
to all orders of perturbation theory. This is familiar from the work of Gawedzki 
and Kupiainen. 
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